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Abstract. - We study the time-dependent height fluctuations of an active membrane con- 
taining energy-dissipating pumps that drive the membrane out of equilibrium. Unlike previous 
investigations based on models that neglect either curvature couplings or random fluctuations 
in pump activities, our formulation explores two new models that take both of these effects 
into account. In the first model, the magnitude of the nonequilibrium forces generated by the 
pumps is allowed to fluctuate temporally. In the second model, the pumps are allowed to switch 
between "on" and "off' states. We compute the mean squared displacement of a membrane 
point for both models, and show that they exhibit distinct dynamical behaviors from previous 
models, and in particular, a superdiffusive regime specifically arising from the shot noise. 



Introduction. - While the physics of biomembranes in equilibrium is fairly developed 
PP, recent studies focus on active membranes that contain proteins, such as ion channels, 
ion pumps, and photo-active proteins like bacteriorhodopsin. These proteins consume the 
chemical energy of ATP, dissipate it into the medium, and thus drive the membrane out of 
equilibrium . The importance of active processes has been demonstrated in an experiment 
shovifing that the fluctuations in the shape of red blood cells depend on the viscosity of the 
environment and on ATP concentrations jjj. In m vitro experiments, nonequilibrium forces 
arising from ion pumps embedded in a membrane are shown to enhance its fluctuations [SUHIE]. 
There are currently two theoretical models for active membranes The Prost-Bruinsma 
(PB) model takes nonequilibrium forces in the form of active noises that include diffusion and 
the stochastic nature (shot noise) of the pumps, but ignores the coupling between the pumps 
and membrane curvature [51 1^ [TU]. The other model proposed by Ramaswamy, Toner and 
Prost (RTP) incorporates this coupling but ignores the random nature of the protein activity 

I12| . For steady state measurements of active membranes, the RTP model agrees quite 
well with experiments In this Letter, we further explore the dynamical properties of the 

RTP model, argue that it is important to include the shot noise for dynamical measurements, 
and present two new models that include both curvature effects and pump stochasticity. In 
the first model, which may be an appropriate description for light-activated pumps such as 
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bacteriorhodopsin, the magnitude of the nonequilibrium force fluctuates on a time scale that 
is fast compared to that of membrane fluctuations. The second model, the two-state model, 
which may be realized in typical ion channels 2 , describes pumps that are able to switch from 
"on" to "off' state on a time scale that is long compared to membrane fluctuation time. These 
models represent natural generalizations of the RTP model, and they show distinct dynamical 
behaviors. In particular, the two-state model predicts that the mean-squared-displacement 
(MSD) of a membrane point exhibits superdiffusion in the experimentally relevant regime, 
whereas the RTP model would predict subdiffusion. Our predictions should be accessible to 
microrheological experiments as carried out in Ref. for similar systems. 

Membrane dynamics in the RTP model. - Consider a 2D tensionless, asymptotically 
flat, fluid membrane embedded a 3D space, lying on the x-y plane. Confined to move in the 
membrane are two types of mobile active pumps. They can either be oriented up or down 
(for the orientation of the force dipole see jj^l) with respect to the normal of the membrane, 
whose shape is described by a height field h{r,t) Their local coarse-grained densities are, 
respectively, p^{r,t) and p^{r,t), where r denotes the 2D position vector. The RTP model 
describes the dynamics of ft,(r, i) and the imbalance field, ip{r,t) = {r,t)—p-^{r,t), taking into 
account the curvature coupling and nonequilibrium forces arising from the pump activities. 
Here, we first summarize the RTP model and further discuss its dynamical properties. 

In the regime where inertial effects are negligible, the dynamics of an active membrane 
is governed by Darcy's law, which states that the relative flow of the solvent with respect 
to the membrane is proportional to the normal force per unit area, Fm{v^t), exerted on it: 
dth{r, t) — Vsz{j^, t) — Api^m(r, t), where Wsz(r, t) is the z-component of the fluid velocity field at 
the surface of the membrane, Xp is the membrane permeability, and Fm{i^, t) is the sum of two 
contributions: a passive part, Fp{r,t) = —STi,/5h{r,t), arising from membrane elasticity and 
an active part from pump activities [El^Ji ^^{1^, t) = Fi[r, t) ip{vt) + F2{v, t) p{r, t) V^h{r, t), 
where p(r, t) = p^(r, t) + p-^{r, t). Note that the expression for Fa{v, t) is the only possible form 
that is linear in the density fields obeying the symmetry: {h,ip, p, Fa) {—h,—^p,p,—FA)- 
In the RTP model, p(r, t) = po is assumed to be uniform, and more importantly, -F'i(r, t) = Fi 
and i^2(r, i) — F2 are assumed to be constant in space and time, capturing only an average 
force produced by each pump. (For simplicity, we set _F2 = throughout this paper, i.e. 
ignoring higher order contributions from the F2 term.) The membrane free energy is given by 



where n is the bare bending modulus, x is the osmotic modulus of the pumps, and S is 
the curvature coupling parameter which arises from head-tail asymmetry of the pumps. In 
equilibrium, k is renormalized to Kg = k — S^/x, which must be positive to ensure stabil- 
ity. The fluid velocity v(x, t) surrounding the membrane obeys Stokes' law: ?7V^v(x, t) = 
Vp(x, t) — Fm (x, t) , where 77 is the solvent viscosity, p{x, t) is the pressure field which ensures 
the incompressibility condition: V • v = 0, and Fm(x, t) is the total force exerted on the fluid 
by the membrane. It contains two parts: a passive part, rp(x, t) = — [5Ti/5h{Y,t)]5{z)i, 
arising from membrane elasticity and an active part arising from the force exerted by the 
pumps on the fluid, modelled as a dipolar force density with force centers located ai z — w 
and z = —w' O E]: Fyi(x, t) — Fa{y, t)[5 {z — w) — 5 [z + w') ] z. Assuming V'(r, t) obeys 
the conserved dynamics: dtip — KV'^ 5Ti / 54> + v and eliminating the fluid velocity from Darcy's 
law, we obtain in Fourier space two coupled Langevin equations 




(1) 



9t/i(q, t) + ujhh{q, t) ^ (3 -!/;(q, t) + /i(q, t), 
5tV'(q, t) + uj^tpiq, t) = 7 h{q, t) -f :/(q, i). 



(2) 
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where tOh = nq^/i-ivi) + nXpq*, (3 = XpFi - {Piw + E)q/{4'q) -^Xpq^, lo^^ Kxq^ = Dq\ D 
is the pump's diffusion constant, 7 — — ASg*, and Pi = Fi{w'^ — w''^)/{2w). The noises 
/x(q, t) and i^(q, i) are assumed to be thermal in origin: {fi{q,t)) — (i^(q, t)) = 0, and 
{^liq,t)^i{~^,t')) = 2kBT[Xp + l/{Ar)q)]S{t ~ t') = TiS{t~t') and {iy{q,t) iy{-q,t')) = 
2kBTh.q^6{t — t') = T25{t — t'), where A is the pump mobihty, ks is the Bohzmann constant 
and T is the temperature. 

It is straightforward to evaluate the two-point correlation function, (ft,(q, <)/i(— q, 0)): 

/i q, t)h{-q, = ^ + ^ ^ + ^ ^ 3 



where A = ujhUJ^4, - f3j, B ^ Uh + cj^, cj± = [ujh + uj^± ^J{uh - t^v)^ + 4/^7 ]/2, and M± = 
Au)± —uj^ ojip. We assume that ijJ± > 0, so that the system is dynamically stable and reaches a 
steady state whose distribution is Gaussian, as expected. Its variance is given by [setting t = 
in Eq. ©], (/i(q, 0)/i(-q, 0)) = [Vi[A + Lol)+ T2p^ ]/{2AB) ^ ksT^s/ {Keq"), in the Ap - 
limit [ig, where T^e = KeT[l + Piw(S + Piw) / {xk)]I {k^ - PiwE/x) HE). This effective 
temperature Tgff is found to be largely consistent with the experimental observations. 

Here, we now discuss the dynamical aspects of Eq. Q. For simplicity, we assume that 
the curvature coupling is small, i.e. uju'^jI) ^ /37, which is indeed the case for the experiments 
in Ref. j7j. However, this approximation may not be always true in general and some of the 
conclusions below may be modified in the strong coupling limit. Within the small curvature 
coupling approximation, Eq. simplifies to 



{h{q,t)H-q,0)) ^ e-^-' 
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2^;, 2(a;2-a;2) 



(4) 



The central quantity of experimental interest is the MSD of a membrane point defined 
by {/^h\t)) ^ {[h{v,t)-h{v,Q)f) = !:llqidq/2n)[{h{q,0)h{-q,0)) - {h{q,t)h{-q,0))], 
where a and L are, respectively, the microscopic and the system size. They introduce two 
time scales ta ~ r}a?/K, and ~ rjL^/n that are assumed to be, respectively, the shortest 
and longest time scales in the problem. Note also that (A/i^(t)) roughly corresponds to trans- 
verse fluctuations of a particle embedded in a membrane in microrheological experiments J7] . 
From Eq. Q), we see that there are two contributions: {Ah'^(t)) = {Ah'^f^{t)) + {Ahl{t)). 
The flrst term comes from thermal fluctuations with the well-known scaling law: {Ahfi^{t)) = 
)t^/3 PHI- The second term arises from the diffusion of the pumps. There 
are two cases to consider; i) in the permeation dominated regime, in which (3 ~ XpFi, we flnd 

r fc^TAA|£fr^ 5/3 , ^ , 

I 4^^1ii(Vici), for 

where tci = K^/{16rfD^). Thus, there is a superdiffusion at short time which crosses over to 
an almost normal diffusion at long time. This is similar to the flndings in Ref. |1(J) . where 
the dynamics of the PB model was analyzed. However, Eq. Q does not contain a term that 
arises from the switching of the pumps that is explicitly included in the PB model. Note also 
that the superdiffusion in Eq. © is a purely nonequilibrium phenomenon in the sense that 
membranes with passive inclusion only show subdiffusion ^19, . ii) In the experimental relevant 
regime in which Xp = 0, i.e. permeation is negligible and (3 = —{Piw + E)/{4:ri) q = ~(3iq, 
we now flnd 

* " \ 1.08 f tor i»f,,. * ' 
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Fig. 1 - The active part of the MSD {Aha(t)) (empty symbols) and the overall MSD (full symbols) as 
a function of r = t/tL for the RTF model in the permeation dominated regime (Ap 7^ 0, Left) and the 
regime with negligible permeation (Ap = 0, Right). The solid lines are scaling laws of Eqs. Q and 
@ which compare well with points obtained from numerically integrating Eq. 0. The parameters 
used are drawn from experimental values in Ref. |H||7|: k = 10 ksT, r) — 10"'^ kg m~^s~^, w = 5nm, 
w' = 4nm, H = wfcsT, Pi = lOksT, D ~ 10"^^m^/s, L = 1mm, po = 3 • lO^'^ m'^, and when 
Ap 7^ we choose /3 = 5 • 10~^*m"^/s. In both cases, the crossover time is td/^L ~ 10~®, where 
tL = riL^/n ~ 10* s. 

In contrast to case (i), the MSD shows an almost normal diffusion at short time and a sub- 
diffuse regime at long times. We have verified these scaling laws [Eqs. (jSJ and ©] by numeri- 
cally computing the MSD as shown in Fig.^ For case (i), the thermal MSD dominates at short 
time, while the active MSD dominates at long time. In constrast, for case (ii), the thermal 
MSD is the dominant contribution. Note also that Eq. @ suggests that membranes contain- 
ing active pumps and those containing passive pumps obey the same scaling law. Therefore, 
the RTF model leads to the conclusion that dynamical measurements cannot distinguish ac- 
tive membranes from passive ones in the experimentally relevant regime. However, as noted 
above, the RTF model ignores the stochastic nature of the pumps, which may be a serious 
approximation made in so far as dynamics is concerned. The natural question to ask is: how 
does the shot noise contribute to the membrane dynamics? We address this question below. 

Active membranes with temporal force fluctuations. - Recall that the active force density 
Fa{y, t) contains Fi (r, t) and F2(r, t) that are assumed to be constant in the RTF model. Thus, 
the most direct way to incorporate pump stochasticity is to allow them to fluctuate in time. 
Here, we focus on the direct force fluctuations: i^i(r,i) — Fi ~\- 5F{t) in the Ap ~ regime. 
Assuming separation of time scales, 5F{t) is described by a Gaussian white noise with zero 
mean: {5F{t)) = and {5F{t)5F{t')) = 2W5{t-t'), where W characterizes the strength of the 
fluctuations. This generalization of the RTF model may capture the effects of fast temporal 
fluctuations of the nonequilibrium forces exerted on the fluids by pumps. Incorporating 5F{t) 
into the RTF equations, Eq. we see that it contributes an additional term i5/3(i) ■0(q, t) 
to the /i-equation with {5P[t)5(3{t')) = Ts5{t - t'), where F^ = 2W[n{q) / {^riq)Y and n{q) = 
(1 + qw) e-«'" - (1 + qw') e"?'"'' is the "structure factor" for the force dipole IHUE]. 

The Langevin equations of this model involve multiplicative and additive noises in a spa- 
tially extended system [201 1^, which in general could pose a mathematically challenging 
problem; however, in this case they can be solved exactly. It is straightforward to derive 
the corresponding Fokker-Flanck equation |22, which describes the time evolution of the 
probability distribution V[{h{q)}, {'ip{q)};t]. Since different q's are decoupled, we can write 
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Fig. 2 - The active part of the MSD {Aha{t)) (open squares) and the overall MSD (solid squares) in 
the presence of temporal force fluctuations as obtained from Eq. For short time, the behavior 
of the active MSD is drastically different from that of the RTF model (open circles). Inset: The 
scale dependent effective temperature Tf,ff{q) induced by force fluctuations (solid curve). In the RTF 



model, T^s/T — 1.9 (straight line) for parameters listed in Fig.^and W = 3.4 ■ 10 



Fig. 3 - Same as Fig. but for the two-state model as obtained from Eq. IIH . The parameters are 
same as Fig.0 and ka = 0.1 s~^, tr ^ 10~^ s, and Fa — 10^* s~^m~^. In contrast to the RTF model 
in the A ~ limit, this MSD exhibits superdiffusion which scales as ~ t^^^. 



P[{/i(q)},{V'(q)};t] ^riq P[Mq),^(q);^], where P[/i(q),V(q);i] satisfies 

dtP = dh [{iohh - 134,) p]+d^ [{uj^ij - jh) P] + (ri/2) dlP +{r2/2) 9^p+(r, v^Vs) d^P. (7) 

The last term in Eq. Q, arising from force fluctuations, is a nonlinear term which may render 
Eq. Q difficult to solve. However, we find, quite remarkably, a closed set of equations for the 
moments of the form ^'m.fe(i) = (/i™(i) 4'^{t)), which directly follows from Eq. Q: 

dt'^m,k{t) = -rn(w;i*m,fc - /3*,„_i,fc+i) - fc (w^*™,fc - 7^'™+i^fc-i) 

+ m(m-l)(ri*„_2,fe + r3*„_2,fc+2)/2 + r2fc(fc-l)^'™,fc-2/2, (8) 

where m = n — fc, 0<fc<n, and n is a positive integer. From Eq. we can compute all 
the moments of the steady-state distribution which, in contrast to the RTF model, is clearly 
non-Gaussian. Furthermore, we can obtain the exact two-point correlation function: 

q, i -q, = ^ ^ + ' ^ 2w T^' 9 

where N± = f3'^uj± + Fs[B {uj± - ojh) + 7/3 ]/2. Note that Eq. Q reduces to Eq. Q if we 
set Fs — 0, as it should be. It is interesting to observe that the denominator in Eq. 
contains the factor AAB — Fs7^, which could become zero at finite q for sufficiently large 
W ^ signaling a dynamical instability. An analysis of this instability will be presented in a 
separate publication, along with the mathematical details leading to Eq. For small W , we 
find that (/i(q, 0)/i(— q, 0)) ~ kBTcf{{q) / {ne<f) , with a scale dependent effective temperature: 
Tcff(g) = KeT[l + Wn{qYl{A7^qx) + Piw{E -f P^w) / {x^)]/ {ne - PiwE/x), which is plotted 
in the inset of Fig. 13 It shows that force fluctuations greatly enhance membrane fluctuations 
mostly in the region q ^ 1/w. Note also that Tcs{q) depends on the viscosity of the solvent, 
whereas Teg in the RTF model does not. Furthermore, temporal force fluctuations contribute 
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significantly to the active MSD, i.e. the second term in Eq. 0, as can be seen in Fig. [3 Note 
that the active MSD in our model is drastically increased at short time, though it does not 
predict superdiffusion. Thus, in contrast to the RTF model, the short time behavior of the 
overall MSD is now dominated by the active contributions. 

Modelling shot-noise with a two-state model. - A two-state model has recently been 
introduced to address the dynamical instability of an active membrane containing inclusions 
with two internal conformational states Such a two-state model of pumps switching 

between "on" and "off' states may also capture the stochastic nature of the pumps. It may 
be a natural model for typical ion channels since they undergo random transitions between 
open and closed states |21 El ■ We model the transition between "on" and "off' as a chemical 

reaction with rate constants kp and ka'. on -5- off on, and introduce imbalance fields for the 
"on" pumps, ipa{r,t) = plnir,t) - p^^{r,t), and "ofF pumps, ipp{r,t) = plg{r,t) - plg{r,t). 
The active contribution to the force exerted on the fluid is assumed to come only from the "on" 
pumps: Fyi(x, i) — Fitpa{i^,t) [S{z — w) — S{z -\- w')] z. In analogy to Eq. the general 
equations of motion for the two-state model in the regime with negligible permeation (Ap = 0) 
are 

dth{q, t) + ujhh{q, t) = /JaV'aCq, t) + Pptppiq, t) + /x(q, i), 

5tV'a(q, t) + WaV'a(q, t) = 7a /i(q, t) + katpp{q, t) + z/a(q, t) + ^(q, t), (10) 
dtippiq, t) + ujptppiq, t) = jp h{q, t) + kpipa{<l, t) + i^p{q, t) - ^(q, t), 

where w/i = Kq^/(Aff), = -(Piw-l-Sa)q/(4?7) = -(iiq, Pp = -Epq/{Ar]), 0Ja,p = Da,pq'^-\-kp^a, 
and 7a, p = — Aa,pSa,pg^. As in the RTF model, ^(q, t), :^a(q, and i'p{q,t) are thermal 
noises with zero mean and {p.{c[,t) q, i')) = ksT / {2-qq)5{t — t') and {vi{(\,t) Uj{—(\,t')) = 
2kBTAi q^Sij S{t—t'). In Eq. H1U|I . ^(r, t) denotes the "chemical noise" associated with the tran- 
sitions between the two states; since the switching process may involve ATF, which is clearly 
not an equilibrium process, its variance is not constrained by the Fluctuation-Dissipation The- 
orem. We simply assume that ^(r, t) has zero mean and {£_{r,t)£_{r',t')) = F3 J(r — r')5(t — t'), 
with F3 being constant in space and time. 

Although an exact, but complicated, expression for the two-point correlation function 
can be obtained analytically, we find it convenient to introduce the following simplifying 
approximations in order to illustrate the essential features: Aa = Ap = A, Xa — Xp — X; 
Da = Dp ~ D, Pa ^ Ppj and LOh^a ^ Pala- Under these approximations, we obtain 



(11) 

where Fi = kBT/{2r]q), T2 = 2kBTAq^, tui = Dq^, 0J2 = Dq^ + 1/tr, and tr = l/(fca + kp). 
The first two terms in Eq. (|ll|l have the same physical origins as the corresponding terms 
in Eq. Q). In particular, the second term gives rise to a MSD which has the same scaling 
laws as in Eq. ©, except that its magnitude is reduced by a factor of ^ (kaTfi)^. The third 
term, which is absent in the RTF model, arises from the noise associated with the switching 
process between "on" and "off' states of the pumps and their diffusion. It is analogous to the 
shot-noise term in the FB model jJUl but ours contains curvature couplings, which are absent 
in the FB model. Note also that the analysis of Ref. TUI assumes that the system is in the 
permeation dominated regime |16|. Assuming F3 is sufficiently large, the MSD arising from 



{h{q,t)h{~q,0)) = 



T2PikaLJp 



2ujh 



PI 
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the third term in Eq. {Ah'^{t)), has to the foUowing scaling laws 



^^'^'W^^ niRrT2 / 4/3 5/3 , (12) 

L 0.16r3/3f(77/K)*/'^T^' , tor t:$>TR. 

Therefore, we find that the MSD exhibits superdifFusion t^^^ at short time. This is qualita- 
tively different from the prediction of RTP model in Eq. (jSJ as well as the force fluctuations 
model in Eq. Since the first two terms in Eq. exhibit only subdiffusion at short time, 
(A/i^(t)} is the dominant contribution to the MSD, which is plotted in Fig. 13 As an estimate, 
TR may be about 0.01 s for typical ion pumps [2 Ell- 
in summary, we have generalized the RTP model by incorporating the shot noise of the 
pumps and demonstrated its significance to the dynamics of an active membrane. 
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